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Description 

[0001] La presente invention concerne une methode pour gen6rer un maillage 3D dans le but de modeliser un corps 
quelconque dont la geometrie est connue et imposee a priori tel que par exemple une structure geologique en tenant 

5 compte fidelement de sa geometrie dans toute sa complexity naturelle. Le maillage realise par la methode s'accorde 
bien dans le cas d'une structure geologique, avec la position des interfaces entre les couches, avec les failles meme 
de forme complexe, les biseaux et autres irregularites de structure, en les representant avec leur conformation reelle. 
[0002] La methode selon Pinventiontrouve des applications notamment pour la moderation en geomecanique, pour 
la caracterisation des reservoirs : couches de sediments, stratigraphie, et aussi pour la simulation des ecoulements 

10 dans les bassins et reservoirs. 

I) ETAT DE LA TECHNIQUE 

[0003] La moderation d'une structure geologique est realisee a parti r d'un ensemble de donnees disponibles ac- 
is quises soit par des mesures, soit par des analyses et des interpretations de donnees sismiques ou des diagraphies 
realisees dans des puits. Ces donnees se presentent sous la forme d'un nuage de points entre lesquels, au mieux, il 
convient d'installer un maillage. La formation d'un maillage en 3D joue un role crucial pour faire des simulations nu- 
meriques dans differentes applications geologiques. Beaucoup de recherches ont ete menees sur le sujet principale- 
ment en ingSnierie de reservoir en raison de ses applications dans le domaine petrolier. 
20 [0004] On connait de nombreuses techniques de maillage 3D pour des simulations numeriques par elements ou 
differences finis et des ensembles de logiciels de simulation crees pour leur mise en oeuvre qui sont adaptes a des 
applications et respectent des contraintes determinees. 

[0005] II est connu pour certaines applications d'utiliser des maillages de type tetraedrique non structure. II est connu 
aussi de faire du maillage hexaedrique de peau ou de surface et de repartir des points dans le maillage soit par 

25 interpolation transfinie, soit encore par resolution ©"equations aux derrvees partielles. Les surfaces sont generalement 
creees a partir de fonctions parametriques et les points du maillage sont obtenus simplement par echantillonnage de 
ces fonctions. Une telle methode ne convient pas generalement pour des applications dans les geosciences ou les 
volumes a modeliser ont des formes le plus souvent tres complexes, diff idles a traduire par des fonctions parametriques 
compte-tenu de la precision qui est recherchee. Aussi prefere t-on approcher les formes complexes par des elements 

30 de surface triangules. 

[0006] Differentes techniques de maillage appliquees a la moderation de reservoir, sont decrites par exemple dans : 

Georges P.L. : "Generation Automatique de Maillages", Editions Masson, 1991 ; 

35 - Guerillot D.R. et al: An ineractive 3D Mesh Builder for Fluid Flow Reservoir Simulation; in SPE - Computer Appli- 
cations, Dec. 1993; 

- Aziz K.; Reservoir Simulation Grids : Opportunities and problems. SPE-JPT, July 1993; ou 

40 - Heineman Z.E. et al : Gridding Techniques in Reservoir Simulation, fntl Forum on Reservoir Simulation. Alpbach 
Austria, Sept. 12-16, 1988. 

[0007] Cependant, la plupart des techniques de maillage connues, et de la les simulations qui en resultent, ne con- 
viennentbien que pour modeliser des structures geologiques relativement simples. Les irregularites geologiques cou- 

45 rantes dans les zones-reservoir model isees, telles que les couches ayant de grandes variations d'epaisseur, les failles 
ou les biseaux, sont rarement pris en compte ou alors seulement sous une forme grandement simplified. Le document 
RENARD P ET AL: 'Three-dimensional geometric modeling of a faulted demain: the Soultz Horst example (Alsace, 
France)' COMPUTERS & GEOSCIENCES, NOV. 1994, UK, vol. 20, no. 9, pages 1379-1390, prend en compte les 
irregularites geologiques en prevoyant un traitement specifique des volumes des couches presentant des irregularites 

so par construction d'une surface triangulee qui passe au mieux par les donnees existantes et dont la geometrie s'accorde 
avec la configuration geologique. 

[0008] Dans les modeles mail les obtenus dans le cadre des geosciences, les failles apparaissent verticalement; au 
mieux, elles sont inclinees et planes. Or il est important de rendre fidelement, la geometrie en 3D des sediments car 
ils possedent des caracteristiques physiques prop res, et aussi la geometrie des failles et autres discontinuites meca- 
55 niques. 

[0009] La methode selon Pinvention permet de g6nerer un maillage 3D respectant la geometrie d'une zone d'un 
milieu complexe dont la geometrie est connue et imposee a priori, tel qu'une structure geologique, dans le but de 
reafiser un modele representatif de ce milieu. La structure est definie a partir de donnees disponibles acquises soit 
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par des mesures, soit par des analyses et des interpretations de donnees sismiques ou des diagraphies realisees 
dans des puits. Elle comporte une modelisation prealable des lirnites (bords de la zone, discontinues ou failles par 
exemple) par des surfaces triangulees, enveloppant au mieux ces lirnites, et un fractionnement de la zone d'etude 
suivant en blocs ayant au plus six faces, dont les faces correspondent a des discontinuites ou des failles principales 
5 de la zone prealablement triangulees. 

[0010] Elle est caracteYisee en ce que'elle comporte : 

la subdivision de chaque bloc en un ensemble de cellules elementaires necessaires pour modeiiser la zone par 
Elements finis, par une reticulation en 3D de chacun des differents blocs obtenus, subdivision qui comprend : 

10 

a) le maillage de chacune des faces du bloc par interpolation a partir de ses bords, de facon a creer une grille 
intermediaire, suivi d'une projection sur la surface trianguiee des noeuds de cette grille interm6diaire suivant 
une direction definie (perpendiculairement a elle-meme par exemple); 

is b) le peuplement en noeuds de grille de chaque bloc dans son epaisseur, par une methode d* interpolation a 

partir des noeuds des differentes grilles formees sur les differentes faces triangulees de chaque bloc, et 

c) un rearrangement dans le positionnement des noeuds des cellules elementaires creees, obtenu (de facon 
iterative, generalement) par relaxation dans le positionnement des noeuds, en fonction de criteres predefinis 
20 tels que la regularite, le choix de la faille relative des differentes cellules et le mode de progression de leurs 

tailles, la direction de projection etc. 

[0011] Chaque bloc etant ainsi reticule en 3D, on peut proceder a un reassemblage ou reempilement de tous les 
blocs. 

25 [001 2] On peut proceder success ivement a une premiere operation de relaxation dans le positionnement des noeuds 
sur chacune des differentes faces triangulees, le peuplement en noeuds internes etant effectue a partir des grilles 
relaxees, et ensuite a une deuxieme operation de relaxation dans le positionnement de ces noeuds internes. 
[0013] De preference, la premiere operation de relaxation est effect uee par des phases iteratives de relaxation- 
projection. 

30 [0014] La methode selon ('invention permet d"obtenir un modele qui respecte trfcs precisement la geometrie souvent - 
complexe de structures geologiques telles que des champs petroliferes, ce qui a pour effet que les simulations qui 
peuvent etre faites a partir de ce modele rendent plus fideiement compte des phenomenes. La methode utilis6e permet 
en outre une economie tres importante de temps de calcul pour realiser les operations de maillage 3D habituellement 
tres longues avec les methodes existantes. 

35 [0015] D'autres caracteristiques et avantages de la methode selon I'invention, apparaTtront a la lecture de la des- 
cription ci-apres de modes de realisation decritsatitre d'exemples non limitatifs, ense reterantaux dessins annexes ou : 

les Fig. 1 et 2 montrent des exemples de relaxation-projection sur un maillage surfacique 4X4; 

40 - la Fig. 3 montre un cube deforme avec des faces triangulees; 

la Fig. 4 montre I'enveloppe externe d'un cube maill6; 

les Fig. 5 et 6 montrent les surfaces internes d'un cube mailie respect ivement avant et apres relaxation; 

45 

la Fig. 7 montre des courbes relevees sur cinq coupes cTune structure geologique de tres grandes dimensions, sur 
laquelle la methode selon I'invention a 6te testee, qui constituent une partie des donnees initiates dont on dispose, 
sur cette structure; 

so - |a Fig.8 montre d'autres donnees initiates sur la meme structure, en Poccurrence, des cartes d'isobathes et les 
emplacements de 4 puits; 

la Fig.9 montre une modelisation surfacique ( base de surfaces triangulees) de la structure geologique 6tudi6e; 
55 - (a Fig.1 0 montre un exemple de decoupage de la structure en blocs; 

la Fig. 11 montre le resultat cfun maillage bloc par bloc de fenveloppe externe de la structure geologique; 
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la Fig. 12 montre le resultat d'un maillage bloc par bloc des surfaces internes de la structure geologique; 

la Fig. 13 illustre a titre cfexemple comment on prend en compte un biseau dans le maillage de la structure geo- 
logique; 

5 

la Fig.1 4 illustre a titre d'exemple comment on prend en compte une faille dans le maillage de la structure geolo- 
gique; 

la Fig. 15 illustre la m&hode cPinterpolation 2D de Coons; 

10 

la Fig.1 6 montre les elements pris en compte dans la methode d" interpolation 3D de Coons; et 

la Fig.1 7 illustre le mode de relaxation-projection des noeuds de maille. 

is [0016] La structure geologique de grandes dimensions illustr6e par les Fig. 7-10. correspond a un dome anticlinal 
comportant six interfaces entre des sediments. 

[0017] La realisation d'un maillage en 3D par elements finis d'une zone a partir d'un ensemble de donnees initiates 
connues pour la zone : interfaces entre diffSrentes couches ou horizons, failles principales, carte en courbes de niveau 
avec les emplacements de puits passe par une succession d'etapes qui seront decrites ci-apres. 

20 

I) Modelisation des limites ou bords de la zone 6tudiee et de ses principales discontinues. 

[0018] On realise une modelisation des surfaces : bords, interfaces, failles en utilisant une logiciel specialise connu 
tel que GOCAD (marque dgposee) sous fa forme d'une juxtaposition continue d'6 laments de surface triangulares. 

25 Dans un premier temps on forme une ebauche qui s'adapte au mieux aux donnees initiales. De maniere iterative, on 
deplace les sommets de la surface en utilisant un logiciel ^interpolation et on densifie localement les elements de 
surface triangulares, de fagon a plaquer au mieux la surface sur le semis de points (donnees d* entree). Les elements 
de surface ainsi crdes ne se raccordent generalement pas comme its le devraient : les interfaces ne sont pas limitees 
aux intersections et au bord du domaine d'etude et elles peuvent se croiser en debordant. On calcule alors les lignes 

30 d' interception entre surfaces et on limite celles-ci aux courbes d' intersection. La Fig 7 montre le resultat final de cette 
modelisation des limites de la structure et f'on peut remarquer que la plupart des triangles des surfaces triangulees du 
biseau sont localises dans les regions de forte courbure. 

[001 9] A ce stade, il n'est pas encore possible de produire un maillage 3D tenant compte des contraintes geologiques 
et numeYiques. Pour y parvenir, on procede a un fractionnement de la zone. 

35 

II) Fractionnement de la zone etudiSe 

[0020] On decoupe cTabord la zone en macro-blocs ayant au plus six faces. II s'agit generalement de blocs de forme 
hexaedrique (a six faces), mais le n ombre de faces peut etre inferieur a six et descendre jusqu'a deux. Dans ce cas, 

40 pour les besoins de la methode, le bloc est quand meme considered comme un block hexaedrique mais ayant des 
faces degenerees (reduites a un point ou a une courbe). Le mode de decoupage est generalement suggere par la 
geometrie des horizons, des principales failles et leur interaction. Pour constituer ces macro-blocs 3D, on a besoin 
parfois d' introduce des frontieres artificielles. C'est necessaire par exemple quand un bloc a plus de six faces natu relies 
et ne satisfait done pas a la regie choisie d'etre au plus hexah6drique. Dans un tel cas, on rSunit deux faces en une 

45 face complexe pour satisfaire le critere, ou bien on introduit des surfaces de coupure artificielles. 

[0021] Un exemple de decoupage de la structure geologique etudiee en macro-blocs est montre a la Fig. 1 0, ou Ton 
a extrapoie la gaille jusqu'a la frontiere de la structure et verticalement projete Pintersection du biseau sur les autres 
surfaces. 

so in) Maillage 

[0022] On realise une reticulation ou maillage en 3D separement de chacun de ces macro-blocs de facon a le sub- 
diviser en blocs §l6mentaires qui vont former les cellules pour la modelisation par Elements finis recherch6es. Cette 
operation est real is 6e en deux temps. 

55 

1) Maillage des faces 

[0023] a) on cree tout d'abord une grille sur chacune des cinq ou six faces de chaque macro-bloc. Pour cela, on 
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forme une grille intermediate par interpolation a parti rdes quatre bords de la face, en utilisant la methode cfinterpolation 
transfinie de Coons dont le principe est rappele ci-apres. 

[0024] Etant donne" quatre courbes param£triques f1(u), f2(u), g1(v), g2(v) (0 2 u, v 2 1) d6finissant quatre bords 
jointifs (Fig. 15), on introduit une formule pour calcuier la surface la plus tendue qui passe par les quatre bords (qui 
5 interpole ces bords) : 

S(u,v) = (l- M )g1(v0 + M52(v) + (l-v)/l(M) + v/2(u)-[(l-<O(l-v)/^0.0) 
10 +0 -«)vP(0.1) + "(l -v)P(1.0) + mvP(U)] * 

[0025] Les sommets P(i J) sont les quatre coins de la surface. 

[0026] On peut alors obtenir un maillage regulier NxM de cette surface par simple echantillonnage en prenant pour 
is sommets les points S(i, j) correspondant aux ui = i/(N-1), vj =j/( M ' 1 )i i=° & N-1 , j=0 a M-1 . La formule 1 restreinte au 
points du maillage devient : 

-[(I -77^)0 - ^(0.0) + 0 - ^ >-j£r WO, ! ) (2) 

25 

[0027] Dans le cas illustnS, les courbes f1.f2.g1, g2, sont obtenues a partir du logiciel GOCAD precitS et sont done 
des lignes brisees. Pour obtenir un maillage regulier MxN, on echantillonne regulierement (en abscisse curviligne ou 
longueur de segment) les courbes f 1 et f2 en N points, et les courbes g1 , g2 en M points, puis on applique la formule 
30 (2). On peut tout aussi bien obtenir un maillage a progression g§om6trique sur chaque direction, simplement en im- 
posant cette contrainte au niveau du reechantillonnage des courbes. 

[0028] Ce maillage interm§diaire sauf cas exceptionnel ne convient pas car il n'y a pas de raison qu'il coincide a. ce 
stade avec la surface qui a prealablement 6te triangulee par le logiciel de surfacage. L'ajustement necessaire consiste 
a projeter les noeuds de la grille intermediate sur la surface triangulee orthogonalement a celle-ci, ou suivant d'autres 

35 directions , de facon a obtenir une nouvelle grille avec des noeuds places sur la surface triangulee. Quand les quatre 
courbes frontieres sont regulierement echantillonnSes, la methode ^interpolation de Coons presente I'avantage de 
produire une grille reguliere en ce sens que sur chaque courbe isoparametrique tous les segments echantillonnes ont 
la meme longueur. Plus g6n6ralement la methode de Coons conserve les rapports de taille entre elements cons6cutifs 
dans une direction donnSes. Ceci permet d'avoir par exemple des maillages reguliers ou respectant certaines pro- 

40 gressions dans la taille des mailles. Cependant la grille obtenue par cette methode connue perd sa regularite (ou le 
critere recherche^ une fois qu'elle est projetee sur la surface triangulee. 

2) Rearrangement des noeuds de face 

45 [0029] Pour retrouverle critere desire dans la distribution des noeuds de la grille projetee, une regularite satisfaisante 
par exemple, il faut deplacer les noeuds du maillage sur la surface triangulee. Pour ce faire, on precede a une relaxation- 
projection de la grille par ajustements successifs comme indiqu6 ci-apres (dans le cas ou le critere est la regularity), 
en relation avec les Fig.1 , 2. 

[0030] Etant donne" un bloc mailte, le but cette relaxation est de rendre son maillage le plus regulier possible au sens 
50 ou n'importe quelle courbe isoparametrique du maillage soit regulierement 6chantillonn6e. Pour ce faire. on d§p!ace 
de facon iterative tout point du maillage de maniere a le positionner au mieux a egale distance de ses deux voisins 
sur chacune des trois directions (u,v,w) de voisinage. Soit P le point qu'on cherche a deplacer. Pour chaque direction 
i, (i = u, v,w) (Fig.1 7), on trouve !e point Pi qui se situe a egale distance des deux voisins de P sur la courbe de direction 
i. (i = u, v.w). Le point P est alors rempIacS par le barycentre des trois points Pi (i = u, v, w). 
55 [0031] Apres avoir d£plac£ tous les points du maillage de la sorte, ceux-ci ne sont plus sur la surface triangulee. On 
les projette alors a nouveau sur la surface triangulee et on n§pete le processus de relaxation -projection jusqu'a atteindre 
la regularity souhaitde. 

[0032] Ceci est illustre" sur I'exemple des Fig. 1 , 2, correspondant a une face du cube deforme" de la Fig. 3 qui est 
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un exemple synthetique penmettant de bien illustrer les differentes etapes de la methode. Cette surface a ete obtenue 
par simple deformation d'une face carree en preservant ses bords ( le bord de la surface reste un carre). Sur ces 
figures, on peut voir trois maillages 4X4. Le premier est un maitlage plat issu de I'interpolation de Coons appliquee 
aux quatre bords (formant un carre). Le deuxieme est le maillage fin dont les noeuds sont obtenus par projection de 
5 ceux du premier maillage. On constate que ce maillage a perdu de sa regularity. Le troisieme est le maillage regulier 
trace en gras et dont les noeuds sont situes sur la surface. II a ete obtenu a partir du deuxieme maillage apres une 
serie de projections et relaxations. La Fig. 4 montre le resultat de ces operations appliquees a I'ensemble des faces 
du cube deform§. 

10 3) Peuplement en noeuds de grille de chaque bloc dans son epaisseur : 

[0033] Pour prolonger en profondeur les grilles formees sur les faces des blocs (Fig.5), on utilise une generalisation 
au volume de la methode d'interpolation transfinie precedente. 

[0034] Etant donne 12 courbes parametriques f1 (u), f2(u ), f3(u), f4(u), g1(v), g2(v), g3(v), g4(v), h1(w), h2{w), h3 
15 (w), h4(w) (0 2 u, v, w 2 1) definissant les 12 bords d'un bloc tel que montre dans la Fig. 16, une generalisation de la 
formule de Coons fournit une parametrisation volumique a I'interieur du bloc. La formule de Coons dans le cas volu- 
mique est la sutvante : 



20 

V(u, v, w) = (! - I ~ w )/!(«) + (l - v)w/2(u)+ wf 3(a) + v(i - w) f4(u) 
-(I -u){\ - w)£l(v)+ ( | - u)wg2{v) + uwg3(v) + it( I - w)g4(v) 
+( l - «X I " If w) + ( I - a)v/i2(iv ) + uvhX w) - v)/i4( v) 
25 +C(«.vmv; 



30 



35 



50 



(3) 



C( u,v , w)=-2[(1 -u)(1 -v)(1 - w) P(0,0,0)+(1 - u)(1 -v) wP(0,0,1 ) 
+(1 -u)v(1 - w) P(0,1 ,0)+(1 - u)v ivp(0,1 ,1 -v)(1 - iv) P(1 ,0,0) 
-v)ivP(1 ,0,1 )+uv(1 - w)P(1 ,1 ,0)+uvivP(1 .1 .1)] 



[0035] Les sommets P(i, j, k) (i, j, k= 0 ou 1 sur la Fig 1 6) sont les 8 coins du bloc. 

[0036] En combinant les equations (1 ) et (3) on aboutit a une expression de la formule de Coons dans le cas volu- 
mique en fonction des 6 surfaces de bord du bloc, S1(u, w), S2(u p w), S3(v,w) t S4(v,w), S5(u t v), S6(u,v). Mais ici, les 
40 surfaces sont aussi supposees etre obtenues par interpolation de Coons a partir de leurs courbes de bord. La formule 
est la suivante : 



V(«, v,w) = -((1 - v)S\(u t v52(u, w) + ( 1 - n)S3(v\ w) + «S4(i\ w) 

+0 - w)S5(u,v) +vuS6f«,v)]+ - C(h,v,h>) 

4 



[0037] Un echantillonnage quelconque de V(u,v,w) fournit alors un maillage EF du bloc. 

[0038] Cette interpolation en profondeur donnerait une distribution reguliere de noeuds si elle etait appliquee a partir 
des noeuds des grilles formes sur les differentes faces du bloc. Ce n'est pas le cas et la distribution obtenue est 
g£ne>alement irreguliere car la methode d'interpolation en profondeur est appliquee aux noeuds report6s par relaxation 
sur les surfaces triangulees, comme decrit dans I'etape 2 precedente. 
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4) Rearrangement des noeuds de bloc 

[0039] Aussi, pour retrouver une r6gularite suffisante (ou respecter un autre critere choisi), la methode selon I'inven- 
tion comporte egalement une relaxation de la grille 3D de la maniere indiquee pr6cedemment. Les Fig .5, 6 montrent 
s le maillage interne du cube dSforme avant et apres relaxation. Le maillage de la Fig. 6, on le note, possede une regularity 
bien meilleure que le maillage de la Fig.5. 

[0040] Quand chacun des blocs precedents de la structure etudiee a et6 convenablement maille en volume, on 
reassemble tous ces macro-blocs apres leur reticulation 3D. 

[0041] La methode qui vient d'etre decrite peut etre utilise>e pour mailler des heterog6neites diverses telles qu'un 
10 chenal ou encore pour mailler radialement un puits devie. 

[0042] On a dScrit des modes de realisation de la methode appliquee au maillage cTune structure geologique. On 
ne sortirait pas du cadre de I'invention, cependant, en utilisant la methode pour le maillage de tout corps ou structure 
dont la geometrie est connue et imposed a priori. 

15 

Revendications 

1 . Methode pour generer un maillage 3D dans une structure dont la geom§trie est connue et imposee a priori telle 
que par exemple une structure geologique definie a partir de donnees disponibles acquises par des mesures, par 

20 des analyses et/ou des interpretations de donnees sismiques ou des diagraphies realisees dans des puits au 

travers de la structure, dans le but de rdaliser un modele repr6sentatif de cette structure, comportant une mode- 
lisation prealable des limites par des surfaces triangulees, enveloppant au mieux ces limites, et un fractionnement 
de la zone en blocs ayant au plus six faces, dont les faces correspondent a des discontinuites ou des failles 
principales de la zone prealablement triangulees, la methode etant caracterisee en ce qu'elle comporte : 

25 

ia subdivision de chaque bloc en un ensemble de cellules eiementaires necessaires pour modeiiser la zone 
par elements finis, par une reticulation ou maillage en 3D de chacun des differents blocs obtenus, subdivision 
qui comprend : 

30 a) le maillage de chacune des faces du bloc par interpolation a partir de ses bords, de facon a creer une 

grille intermediate, suivi cTune projection sur la surface triangulee des noeuds de cette grille intermediaire 
suivant une direction definie, de facon a g6n6rer des grilles de face; 

b) le peuplement de chaque bloc en noeuds internes par interpolation a partir des noeuds des differentes 
35 grilles formees sur les differentes faces triangulees de chaque bloc, de facon a obtenir une distribution 

r6guliere des noeuds dans le volume de chaque bloc; et 

c) un rearrangement dans le positionnement des noeuds des cellules eiementaires creees, obtenu par 
relaxation dans le positionnement des noeuds, en fonction de criteres predef inis. 

40 

2. Methode seton la revendication 1 , caracterisee en ce que Ton procede successivement a une premiere operation 
de relaxation dans le positionnement des noeuds sur chacune des differentes faces triangulees, le peuplement 
en noeuds internes etant effectue a partir des grilles relaxees, et ensuite une deuxieme operation de relaxation 
dans le positionnement de ces noeuds internes. 

45 

3. Methode selon la revendication precedente, caracterisee en ce que Ton realise la premiere operation de relaxation 
par des phases iteratives de relaxation-projection. 

4. Methode selon Tune des revendications 1 ou 2 ou 3, caracterisee en ce Ton choisit comme direction de projection, 
so une direction perpend iculai re a la grille intermediaire. 

5. Methode selon Tune des revendications pr^cedentes, caracterisee en ce que Con utilise un mode cf interpolation 
de type transfini. 

55 6. Methode selon Pune des revendications pr6cedentes, caracterisee en ce que Ton rearrange le positionnement 
des noeuds des cellules eiementaires creees de facon a obtenir des maillages reguliers. 

7. Methode selon Tune des revendications 1 a 5, caracterisee en ce que Ton rearrange le positionnement des noeuds 
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des cellules 6lementaires crepes de facon a obtenir un mode determine" de progression de la taille des mailles. 

8. Methode selon Tune des revendications precSdentes, caracterisee en ce que Ton procede a un r6assemblage 
de tous les blocs. 



Patentansp ruche 

1. Verfahren zum Erzeugen eines 3D-Netzwerkes in einer Struktur, deren Geometrie bekannt und a priori vorge- 
10 schrieben ist, wie zum Beispiel eine geologische Struktur, die aus verfugbaren Daten definiert ist, welche durch 

Messungen, durch Analysen und/oder Interpretationen seismischer Daten Oder Diagraphien aufgenommen sind, 
welche in Bohrlochern entlang der Struktur realisiert sind, mit dem Ziel, ein fur diese Struktur representatives 
Modell zu verwirklichen, umfassend eine vorherige Modeliierung der Grenzen durch aus Dreiecken zusammen- 
gesetzte Oberflachen, welche diese Grenzen am besten umhullen, und eine Fraktionierung der Zone in Blocke 
15 mit hochstens sechs Flachen, deren Flachen Hauptdiskontinuitaten Oder Hauptverwerfungen der zuvor aus Drei- 

ecken zusammengesetzten Zone entsprechen, wobei das Verfahren dadurch gekennzeichnet ist, dass es um- 
fasst: 

die Unterteilung jedes Blockes in eine Zusammenstellung von Elementarzellen, die notwendig sind, urn die 
20 Zone durch Endelemente durch eine Vernetzung oder 3D-Netzwerk von jedem der verschiedenen erhaltenen 

Blocke zu modellieren, welche Unterteilung umfasst: 

a) Das Netzwerk von jeder der Flachen des Blocks durch Interpolation ausgehend von seinen Grenzen 
derart, dass ein Zwischengitter erzeugt wird, gefolgt von einer Projektion der Knoten dieses Zwischengit- 

25 ters auf die aus Dreiecken zusammengesetzten Flache gemaS einer definierten Richtung, derart, dass 

Flachengitter erzeugt werden; 

b) die Population jedes Blockes mit inneren Knoten durch Interpolation aus Knoten von verschiedenen 
gebildeten Gittern auf die verschiedenen aus Dreiecken zusammengesetzten Flachen jedes Blocks, der- 

30 art, dass eine regelmassige Verteilung der Knoten in dem Volumen jedes Blockes erhalten wird; und 

c) Einer Wiederanordnung in derPositionierung der Knoten dererzeugten Elementarzellen, erhalten durch 
Entspannung in der Position ierung der Knoten in Abhangigkeit vorgegebener Kriterien. 

35 2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dass man aufeinanderfolgend zu einem ersten Entspan- 
nungsvorgang in der Positionierung der Knoten auf jeder der verschiedenen aus Dreiecken zusammengesetzten 
Flachen ubergeht, wobei die Population mit inneren Knoten ausgehend von den entspannten Gittern durchgefuhrt 
wird und man anschlieBend zu einem zweiten Entspannungsvorgang in der Positionierung dieser inneren Knoten 
ubergeht. 



40 



3. Verfahren nach dem vorhergehenden Anspruch, dadurch gekennzeichnet, dass man den ersten Entspannungs- 
vorgang durch iterative Relaxations- Projektions-Phasen durchfuhrt. 

4. Verfahren nach einem der Anspruche 1 oder 2 Oder 3, dadurch gekennzeichnet, dass man als Projektions- 
Richtung eine Richtung wahlt, die senkrecht zum Zwischengitter liegt. 

5. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man einen Interpolati- 
onsmodus vom transfiniten Typ verwendet. 

6. Verfahren nach einem der vorhergehenden Ansprfiche, dadurch gekennzeichnet, dass man die Positionierung 
der Knoten der erzeugten Elementarzellen derart wiederanordnet, dass regelmaBige Netzwerke erhalten werden. 

7. Verfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet, dass man die Positionierung der Knoten 
der erzeugten Elementarzellen derart wiederanordnet, dass ein vorgegebener Progress ionsmod us des Umfangs 
der Maschen erhalten wird. 

8. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man zu einem emeuten 
Zusammensetzen aller Blocke ubergeht. 
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Claims 

1. Method for generating a 3D mesh in a structure whose geometry is known and imposed a priori, such as for 
example a geological structure defined from available data acquired by measurements, by analyses and/or inter- 
pretations of seismic data or logs produced in bore holes through the structure, with the object of producing a 
representative model of this structure, comprising prior modelling of the limits by triangulated surfaces, enveloping 
these limits as well as possible, and splitting of the zone into blocks having at most six faces, the faces of which 
correspond to previously triangulated major faults or discontinuities in the zone, the method being characterised 
in that it comprises: 

subdivision of each block into a set of elementary cells necessary to model the zone by finite elements, by 
creating a 3D lattice or mesh from each of the different blocks obtained, said subdivision comprising: 

a) creating a mesh from each of the faces of the block by interpolation from its edges so as to create an 
'5 intermediate grid, followed by projection onto the triangulated surface of the nodes of this intermediate 

grid in a defined direction so as to generate face grids; 

b) population of each block with internal nodes by interpolation from the nodes of the different grids formed 
on the different triangulated faces of each block so as to obtain a uniform distribution of the nodes in the 

20 volume of each block; and 

c) rearrangement in the positioning of the nodes of the elementary cells created obtained by relaxation in 
the positioning of the nodes as a function of predefined criteria. 

25 2. Method according to claim 1 , characterised by successively a first operation for relaxation in the positioning of 
the nodes on each of the different triangulated faces, the population with internal nodes being effected from the 
relaxed grids, and then a second operation for relaxation in the positioning of these internal nodes. 

3. Method according to the preceding claim, characterised in that the first relaxation operation is carried out by 
30 iterative phases of relaxation and projection. 

4. Method according to one of claims 1 or 2 or 3, characterised in that a direction perpendicular to the intermediate 
grid is chosen as the direction for projection. 

35 5. Method according to one of the preceding claims, characterised in that a mode of interpolation of transfinite type 
is used. 

6. Method according to one of the preceding claims, characterised in that the positioning of the nodes of the ele- 
mentary cells created is rearranged so as to obtain uniform meshes. 
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7. Method according to one of claims 1 to 5, characterised in that the positioning of the nodes of the elementary 
cells created is rearranged so as to obtain a given mode of progression for the size of the mesh. 

8. Method according to one of the preceding claims, characterised in thataW the blocks are reassembled. 
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Fig. 1 
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Fig. 2 
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Fig- 5 
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Fig. 6 
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Fig. 8 
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FIG.16 




FIG.17 
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